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Abstract 

We calculate the leading twist valence quark distribution in the pion in 
the framework of QCD sum rules with nonlocal condensates. Particular 
attention has been paid to the correct account for the bilocal power 
corrections. 

The determination of parton distributions is, up to now, reserved to experimental 
studies but as a final goal they are expected to be evaluated from the first principles 
of the theory. In the lack of complete understanding of the yet unclear confinement 
mechanism they provide a challenging task for nonperturbative methods presently 
available. Among the approaches which account for nonperturbative effects the most 
close to QCD perturbation theory are the QCD sum rules [1]. In the last decade 
they were applied with moderate success to determine nucleon and photon structure 
functions in the region of intermediate values of the Bjorken variable [2] and, recently, 
in the small A-region for the light-cone position representation [3]. While the nucleon 
structure functions are now well defined by the analyses of the precise experimental 
data and are attacked theoretically, much less is known about the parton distribution 
of other hadrons, in particular of vr-meson. Being of interest in their own right they 
provide good testing ground for predictions of the QCD sum rule method. 

In order to evaluate the quark distribution in the pion by means of the QCD sum 
rules method, we consider an appropriate three-point correlation function 

W,.{p„P2,q) = i\j d'xd'ye'-'^^^+'^y{0\T |j^(x), (9 (^y + ^n, y - (0)| |0). 
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of two axial currents that have non-zero projection onto the pion state being proportional 
to the pion decay constant (0|j|^ |7r(p)) = ifjrP^ and the nonlocal string operator O on 
the light cone defined by equation 
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and $ is a path ordered exponential in the fundamental representation of the colour 
group along the straight line C which insures the gauge invariance of the parton 
distribution. This operator when sandwiched between the pion states defines the 
twist-2 valence quark distribution in a hadron [4] 

{h{p)\ol^n,-^n] \h{p))=i f\os{Xx)uH{x). (3) 







The usual strategy is to use the duahty between the hadronic and partonic repre- 
sentation for the correlator under investigation. 

On the one hand, we should consider the dispersion relation for the latter and 
extract the contribution due to the low lying hadron, namely, due to vr-meson, 
approximating the higher state contribution by perturbative spectral density. We 
perform the double Borel transformation and put the parameters equal = M| = 
2M^ in order not to introduce the asymmetry between the initial and final pion states 
and to make contact with two-point sum rules for the pion decay constant. 

On the other hand, we consider the OPE for the same quantity. Of course, the 
QCD sum rules with local condensates are inappropriate here because the usual local 
power corrections produce 5-type contribution to the distribution function. It is 
not surprising since some propagators are substituted by constant factors that do 
not allow the momentum to fiow and the whole hadron momentum be carried by 
a single quark. The probability density of this configuration in the phase space is 
5(1 — x). Higher condensates produce even more singular terms. However, this 
singular contribution can be smeared over the whole region of the momentum fraction 
from zero to unity by avoiding the Taylor expansion of the generic nonlocal objects 
which are the starting point of all QCD sum rule calculations and introducing the 
concept of nonlocal condensate which assumes the finite correlation length for the 
vacuum fiuctuations. 

At the two-loop level, to which we restrict our analysis, we need the bilocal quark 
and gluon condensates, trilocal quark-gluon condensates and four-quark condensates. 
The latter will be factorized into the product of bilocal scalar quark condensates 
via the vacuum dominance hypothesis. For explicit calculations, it is convenient to 
parametrize the bilocal condensates in the form of the well-known a-representation 
for propagators [5] 

(0|V'(0)$[0,x]V'(x)|0) = (V'V') / dafsia)e"'^ l\ 
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We use the following ansatz for the distribution of vacuum quarks in the virtuality a 

Here = 0.2GeV^^ and 7 is fixed from the lowest nontrivial moment of the distribution 
function fs that is related to the value of the average virtuality of the vacuum quarks 

Conventional calculations of the perturbative diagram with the light quark masses 
neglected result in 

^pUM\ Q\ A) = -^M^ dx co^{\x)x-x exp f---^") • (6) 



Note that we have kept the t-channel momentum transferred to be nonzero. If we 
expand the cosine in the Taylor series and integrate over x, we find out that each 
moment possesses logarithmic non-analyticities of the type {Q^)^lnQ^. These terms 
come from the small-x region, where the spectator quark carries almost the whole 
momentum of the pion, so that the struck quark becomes wee and can propagate 
over large distances in the t-channel. Therefore, we have to perform additional 
factorization for separation of small and large distances in the corresponding invariant 
amplitude; this will lead to the appearance of additional terms in the OPE for the 
three-point correlation function which correct the small-x dependence of the parton 
density. 

The simplest nonperturbative correction comes from the vector condensate 
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The dominant contribution is due to the four-quark condensate. Performing 
straightforward calculations we obtain 
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The gluon as well as trilocal quark-gluon condensate contributions are numerically 
much less important than the power correction we accounted for; therefore, we neglect 
them in what follows. 

It is well known that there exists a parton sum rule that implies that the pion 
contains one u-quark. Summing the calculated contributions and taking the formal 
limit — > in the perturbative term we can convince ourselves, comparing the result 
with the sum rule for the pion decay constant [5], that the normalization condition 
is broken. The reason for this has already been mentioned earlier and we elaborate 
this point below. 

As we have seen, in the limit — > the perturbative term though finite contains 
the logarithmic non-analyticities at this point. This is a typical example of the mass 
singularities in the QCD sum rules framework [7, 8, 9]. In order to get rid of this 
perturbative behaviour and replace it by a physical one, it is necessary to modify the 
original OPE. Eor the form factor type problem a two-fold structure of the modified 
OPE has been realized in refs. [10, 11] being of the following schematic form: 

Wiplplq') = j:c^'\plpll'){Od) + Y: I ^/W^-cW(x)WXg,x,A). (9) 



An additional second term determines the contribution due to the long-distance 
propagation of quarks in the t-channel. Here W* are the two-point correlators of 



the operator in question and some nonlocal string operator of a definite twist [9] 
which arises from the OPE of T-product of pion interpolating fields. The coefficients 
C'^'^^pI^pI^ q^) in eq. (9) are free from non-analyticities or singularities in because 
they are defined as the difference between the original diagram and its factorized 
expression which is the perturbative analogue of the corresponding bilocal correlator. 

The simplest bilocal power correction is given by the convolution of the quark 
propagator C^^^(x) = 25*+ (x) and 

X, An) = z Jd'ye'^y{0\T |u(0)7+$[0, x]u(x), O (^y + ^n,y - ^n^ | |0). (f 0) 

We extract the contact term due to the vector condensate from this correlator and 
saturate the remaining part by the contributions of the mesons of increasing spin; 
these are /9°, g states and so on. The net result for the difference between the "exact" 
bilocal and its perturbative part reads 

¥sUm\\) = —iraJuuf dxcos(\x)xfv(xM^) + Q^e^M2 V (iXf-^ (11) 
81 -^0 J=l,3,... 



where Zj = ctj/AM^^ ctj is the continuum threshold and niM is the mass of the 

(2) 

lowest meson state in the channel of given spin J. Here (pj are the wave functions 
describing the light-cone momentum fraction distribution of quarks inside mesons 
that parametrize (to the leading twist accuracy) the matrix elements 

(0|V^(0)$[0,x]7,V(x)|M^(g,ry)) = ea,,...,,_,x,,x,, ...x,,_, {mMY 

where J is a spin of the meson, rj its polarization and e[i'Jj\^2 a^j i ^ polarization 
tensor. 

The former term in eq. (12) is a contact-type contribution due to the vector 
condensate. The first one in the curly brackets is the difference between the perturbative 
analogue of the bilocal correlator and the continuum contribution into the "exact" 
one. This part cancels the logarithmic non-analyticities in the perturbative diagram 
(eq. (6)) corresponding to the leading twist-2 operator in the OPE of pion currents. 
The tower of the next-to-leading non-analyticities can be subtracted in a similar 
way by accounting for the twist-4 operator. The last term displays the physical 
contribution to the correlation function that possesses the correct behaviour in the 
"momentum transferred Q^"- 

The sum of eqs. (6) and (12) is an analytical function in as all singularities are 
replaced by the combination -\- ctj which is safe in the limit — > 0. Due to the 



presence of the non-analyticities in each moment of the distribution function, we need 
an infinite number of parameters to be found from additional sum rules. Obviously, 
this is an impossible task. Safely, for the problem at hand, this part vanishes in the 
forward limit and the sum rule is dominated by the contact terms. 

The dominant contribution comes from the bilocal correlator convoluted with a 
three-propagator coefficient function. All nonperturbative information is accounted 
for in the correlator 



W'liq, X, An) = I d''ye"i^O\T I u(0)$[0, x]u(x), y + -n, y - -n ) }> |0). (13) 
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Extracting the contact term and performing the integration by using the method 
outlined at the beginning of the paper, we obtain the following contribution to the 
structure function: 



^^]^(M^,A) = — 7ras(uu)^ / dxcos{Xx 






Now, having accounted for additional terms in OPE, we can easily check that the 
normalization condition for the quark distribution in the pion is restored. 

Eor zero the perturbative spectral density is concentrated on the line Si = S2, 
so that there is no transition between the states with different masses. We collect all 
contributions and make the continuum subtraction that results in the substitution 
— > M^(l — exp[— Sq / M'^)) in the perturbative term. We have found good stability 
of the distribution function with respect to the variation of the Borel parameter 
in the region 0.5 < < 0.8 for the standard value of the continuum threshold 
So = 0.7GeV^. The normalization point of the OPE is /i^ ~ O.bGeV^; therefore, 
the function obtained can be regarded as an "input" quark distribution at this low 
energy scale. In Eig. 1, we present the curves for the valence quark distribution 
in the pion for = O.GGeV^: the solid and long-dashed lines correspond to the 
values of the average virtuality of vacuum quarks A^ = O.GGeV^ = O.lbAGeV^) 

and A^ = OAGeV^ = O.OSTGeV^), respectively. In the large-x region, the 

corrections due to the quark condensate do not exceed 30% of the perturbative term. 
However, in the small-x region at x = 0.2 the ratio of the contact term to the 
main one comprises 50% for A^ = O.GGeV^ and 70% for A^ = OAGeV^. Below 
this point the nonperturbative contribution increases and reaches 100% at x = 0.13 
for A2 = 0.4^6^2 (for A^ = O.QGeV^ it still amounts 50%). So, for x as small as 0.2 
we could not trust x-dependence of our result. Of course, there is no possibility to 
reproduce the correct x — > behaviour of the parton density in the present approach 
as it is determined by the exchanges of the Regge trajectories. 

Since our result is valid only in the limited region of Bjorken variable, we could 
not evolve it to the experimentally accessible energies. In Eig. 1, we compare 



our calculation with the distribution obtained in the NJL model [12] at the same 
normalization point and find reasonable agreement between two approaches in a wide 
region of the momentum fraction. In Fig. 2, the latter evolved up to = 20GeV^ 
(short-dashed curve) is compared with the presently available fits of experimental 
data. It shows good agreement with the result of the analysis of Sutton, Martin, 
Roberts and Stirling (solid curve) [13], which is consistent with all present Drell-Yan 
and prompt photon ttN data. We also present the (long-dashed) curve due to Gliick, 
Reya and Vogt [14]; however, their result does not agree with E615 experiment which 
requires the valence distribution to be larger by 20%. If the GRV curve is renormalized 
within a factor of 1.2 — 1.3 in the central region, there will be no disagreement between 
the different analyses. 




Fig. 1: Quark distribution in the pion at the Fig. 2: The experimental fits of the valence 

low energy scale /^^ ~ O.bGeV^ calculated M-quark distribution in the pion at = 

from the QCD sum rule for different values of 20GeV^: SMRS result [13] is depicted by 

the average virtuality of vacuum quarks: solid solid curve, GRV analysis [14] is shown by 

and long-dashed curves correspond to = long-dashed one. Short-dashed curve is the 

0.6GeV^ and = OAGeV^, respectively. M-quark density calculated in the NJL model 

Short-dashed curve is the M-quark density [12] evolved up to = 20GeV^. 
found in the NJL model [12]. 

In conclusion, we have calculated the pionic parton density at low momentum 
scale in QCD sum rules with nonlocal condensates. It is shown that the parton sum 
rule is fulfilled only after the bilocal power corrections are accounted for. We have 
found good agreement with the u-quark distribution function computed in the NJL 
model which when evolved up to the experimental scales is well comparable with 
data. 
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